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I. INTRODUCTION 



The meson light-cone distribution amplitude (DA) provides essential information on the 
nonperturbative structure of mesons. In particular, the leading-twist meson DAs play a 
role of input for describing hard exclusive reactions due to factorization theorems 0, 0, 
H 0|. In particular, the pion DA has been investigated extensively in various theoretical 
approaches: For example, in the QCD sum rules (QCDSR) 0,11 00, Hi, in lattice QCD 
(LQCD) 10, , in the chiral quark model (xQM) from the instanton vacuum 

EBB, 

in the NJL models H H 12] an d so on. Recently, the CLEO experiment |l8[ has measured 
the 7*7 — > 7r transition form factor, which gives criteria for judging the existing theoretical 
pion DAs 1^ . Bakulev et al. 2(1 21] have carried out similar analyses. 



In the previous work |1J], we have investigated the leading-twist pion and kaon DAs 
within the framework of the nonlocal xQM from the instanton vacuum with SU(3) sym- 
metry breaking explicitly taken into account. However, it is well known that the nonlocal 
interaction breaks the conservation of the vector and axial- vector currents jl^, HI 0]. 
The nonlocal xQM from the instanton vacuum suffers from the same difficulty due to the 
zero-mode approximation |24[. Since the leading- twist meson DAs involve the axial- vector 
operator, one has to deal with the problem of the current conservation. While Ref. 24] 



proposed a systematic way to remedy this problem of the current conservation, one has to 
handle the integral equation. Refs. pa . Hif derived the light-quark partition function in the 
presence of the external gauge fields. With this gauged partition function, it was shown 
that the low-energy theorem for the transition from two-photon state to the vacuum via the 
axial anomaly was satisfied [25j]. Moreover, the magnetic susceptibility of the QCD vacuum 
was properly obtained ji^ . 

The conservation of the axial- vector current concerning the leading-twist pion DA was al- 
ready discussed in Refs. jl?], by using the instanton-motivated separable nonlocal kernel 
for the effective action and the NJL model with the Pauli-Villars regularization, respectively. 
It was pointed out that the nonlocal contributions are about 30% of the local one j27|. In 
the present work, we want to investigate the leading-twist pion and kaon DAs, employing 
the method developed in Refs. [H |26| with the conservation of the axial-vector current 
considered. We will show that the nonlocal part of the pion DAs arising from the current 
conservation plays an important role in improving the results of the corresponding Gegen- 
bauer moments. 

The present work is organized as follows: In Section II we show how to make the nonlocal 
effective chiral action gauge-invariant from the instanton vacuum. In Section III we show 
how to calculate the leading-twist pion and kaon distribution amplitudes in the present 
approach. In Section IV we discuss numerical results. In Section V we summarize the 
present work and draw conclusions. 



II. GAUGE-INVARIANT EFFECTIVE CHIRAL ACTION 



In this Section, we briefly review how one can make the low-energy effective QCD partition 
function from the instanton vacuum in the presence of external axial- vector fields, following 
the approach developed in Ref. 25|, 26||. The partition function derived from the instanton 
vacuum [H |2t| can be written in Euclidean space: 

exp / d 4 x \if)f(x)(ifi + irrif)il)f(x) 
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+l 1 7^ e ^M(*)>/^ u SV^W^(p)]' (!) 

where m/, and denote the quark fields, the current quark mass, and the momentum- 
dependent dynamical quark mass with flavor /, respectively. The U 15 or M. represents the 
pseudo-Goldstone field defined as follows: 



^ = U{x)^ + U\x)^ = 1 + — l5 M a X a - ^M 2 ■ • • , (2) 



where A stands for the 3x3 Gell-Mann matrix. Musakhanov , Ey| improved the partition 
function in Eq. ([TJ by taking into account effects of flavor SU(3) symmetry breaking, so that 
the dynamical quark mass acquires the contribution of the rrif corrections: 



M f (k) = M F 2 (k) 



d 2 d 



(3) 



where M is the dynamical quark mass with zero momentum transfer in the chiral limit. Its 
value is determined by the saddle-point equation: Mq — 350 MeV. F(k) is the momentum- 
dependent part which arises from the Fourier transform of the fermionic zero-mode solutions 
in the instantons. Note that we parameterize Mo(k) as a simple-pole type to make the 
analytic continuation to Minkowski space easy. 

The momentum-dependent dynamical quark mass Mf(p) breaks the conservation of the 
Nother currents, i.e. it violates Ward-Takahashi identities. We want to show how to remedy 
this problem in the following. We first define the total quark propagator S in the presence 
of the instanton ensemble A and external axial- vector field a M = a^A a /2 and the quark 

propagator S{ with a single instanton A4 as well as the a^. 

§ - ? ^ = - , (4) 

ip + gA + 4l5 + im % ip + gA { + ^75 + irh ' 

where m is the current quark mass matrix: diag(m u , mj, m s ). We assume that the total 
instanton field A may be approximated as a sum of the single instanton fields, A = YliLi ^-ii 
which is justified with the average size of instantons p m 1/3 fm and average inter- instanton 
distance R ~ lfm [i^, H^]. Defining the quark propagator S with external fields a M and 
the free one So as follows: 

^° = -a , 1 1 ~> ^0 = — , (5) 

ip + (2175 + im ip + vm 

the S can be expanded with respect to a single instanton: 

3 = S + J2(3i - So) + J2(3i ~ So^oHSj - S ) + ■ ■ ■ . (6) 

In order to specify the gauge dependence, we rewrite Si and So in the following form: 
Sj = LjS \L ■ . S? 



' ' ' ' ip + gAi + ^ + tm'' 

1 



So - US' LT\ S' 0i - + ^ + (7) 
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where 

4'-y B = L-\ip + ^)L. (8) 
The L denotes the Wilson gauge connection which is expressed as the path-ordered exponent: 

Li(x,Zi) =Pexp^i jT d^a^Ts) , (9) 

where Zj denotes an instanton position. Having carried out the manipulation as in Refs. ^2^, 
Efil ]. we arrive at the low-frequency part of the fermionic determinant: 

Det low = (det(ip + ^ + im))- 1 f J[ D^ f D^\ (10) 

/ 



x e 



where 

V±,f[^fM = J d 4 x(^ f (x)L(x,z)i0 ±fi (x;£ ± )) J d 4 y (<t>l fi (y;( ± )(ipL(y, z)^ f (y)) (11) 

with the fermionic zero- mode solutions $±,0- 

The gauge connection L has some arbitrariness due to its path dependence. However, 
we can show that such dependence can be minimized. To be more specific, we consider the 
extended zero mode: 

(ip + /i + t<y s )\$ o ) = 0, l^) ~ |$ > - W7 5 |$o>, (12) 

where A denotes an instanton field located at z, l^^) stands for the solution in the presence 
of the external field, and S^z the well-known non-zero mode of the propagator in the 
instanton field (see the review [32| and references therein). Here, we have assumed tacitly 
that the momentum of the external field is small. Inserting the gauge connection L(x, z) 
into Eq. ()12j) . we obtain: 

(^ + A + ^75)|$o> = 0, a'^ = L-\id, + a^)L, \%) = L~\x, z)\$ ). (13) 

If we utilize the relations (ift + A)\§ Q ) = and S NZ (ift + A) = 1 - |$ )($ |, the solution 

l^g 1 " 1 ) in Eq. (|T3|) can be reduced to the corresponding solution 1$^) without any problem 
of the path dependence arising from the gauge connection L. 

However, if the zero-mode approximation Snz ~ \ = Soo is used to find the solution 1 

if 

l$oo'> = l$o> - W7 5 |$o>, (14) 
then 

l$ffi> = l$o> - -WtsI^o) + 5 00 {% J ^WOtb) \%) (15) 



1 In fact, the zero- mode approximation is the origin of the problem of the current non-conservation from 
the instanton vacuum. Pobylitsa suggested a more reliable approximation to restore the current conser- 
vation 0. 
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with applying the inverse gauge connection L 1 to Eq. We can easily see that Eq. ([13)1 
depends on the path of the gauge connection. We need to choose the path in such a way 



that || - $^|| 2 can be minimized. As shown explicitly in Ref. |26j, the straight line 
provides the most optimized path. 

Having treated the gauge problem in the presence of the external axial-vector field, we 
finally obtain the low-energy partition function of the nonlocal xQM from the instanton 
vacuum: 



i(i)ii2 



(16) 



Z[a,m] = exp / d 4 x\ijjj(x)(ift + ^75 + imf)ipf( 

d A kd A p 



+ % 



-i(k-p)-x i f , 



with the effective chiral action in the presence of the external axial-vector source 
S e g[a,m] = — iV c Trln 



ip + fas +im f + iyj M f {id^ + a^ 5 ) JJj g M g (id^, + a^js) 



(17) 



where Tr denotes the functional trace and traces over flavor and spin spaces, generically. 



III. MESON DISTRIBUTION AMPLITUDES 

The leading-twist pseudoscalar meson light-cone distribution amplitudes (DAs) can be 
defined with the conserved axial-vector current as follows (here, we consider the light-cone 
gauge A ■ n = 0): 

(0\<lA™)j^ 5 q a (-Th)\M(P)) = iV2F M P, f 1 due^ 2u -V p ™Mu), (18) 

Jo 

where M. denotes the flavor SU(3) octet pseudoscalar meson field with the on-mass shell 
momentum P 2 = m 2 M in the light-cone frame. In the present work, we choose J\A = 
7r + or JW = K + . u and n stand for the longitudinal momentum fraction and light-like 
vector, respectively. Fm is the pseudoscalar meson decay constant. Note that, when the r 
approaches to zero, Eq. (fTH|) becomes the pseudoscalar meson decay amplitude due to the 
DA normalization condition: 
l 

du4>M(u) = 1. (19) 



o 

This equation will be used for determining the scale parameter of the model. At this 
point, we want emphasize that in order to calculate the light-cone DA, it is necessary for 
us to work in Minkowski space in spite of using the instanton framework which is well- 
defined in Euclidean space as discussed before. For this purpose, we assume simple anal ytic 
continuation (the Wick-rotation) between the two spaces in a practical point of view [12|, 



The expression for the leading-twist pion and kaon DAs in the present approach is ob- 
tained as follows: 



Nrh^ r d 4 k 



M f (k) jM g (k-P) 



D f (k) D f (k-P) 
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+ 



M f (k)jM g (k - P) JM f (k)JM g (k - P) 



D f {k) 



Df(k-P) 



(20) 



where JMf(k)p = dJ Mf{k)/dk il . tr 7 denotes the trace over the Dirac spin space. Note that 
the last two terms that we call the nonlocal contributions in the bracket make the present 
result gauge-invariant. The quark propagator Df(k) is expressed in terms of the momentum- 
dependent dynamical quark mass Mf(k) and current quark mass mj as follows |3lj : 



D f (k) =ft-[m f + Mf(k)]. 

We employ a simple-pole type parameterization for F(k) in Eq. (j3J): 

nA 2 



F(k) 



(nA 2 - k 2 + 



(21) 



(22) 



since it is easy for the analytic continuation to Minkowski space. 

The derivative of the dynamical quark mass a/M m appearing in the nonlocal contributions 
can be evaluated analytically with the simple-pole type form factor of Eq. 



dM f {k) _ 2n v /M / (fc) 



2JM f {k) dk„ (nA 2 - k 2 ) 



(23) 



Inserting Eqs. (|2*2*|l and (J2*3|l into Eq. (|2*0*jl . we finally obtain the expression for the leading- 
twist light-cone pion and kaon DAs: 



''2F 2 M 



d A k 
2nk. 



S[uP^ 



k + ]y/M f (k)y/M f (k-P) te, 
2n(k + - P+) 



1 



D f (k)(nA 2 - k 2 ) D g (k - P)[nA 2 - (k - P) 2 } 

t [t l + r NL ] , 



[D f (k) ! '°D g (k-P) 



75 



2^ 



(27TY 



(24) 



where p + stands for p ■ n. The subscripts L and NL stand for the local and nonlocal 
contributions, respectively. Further evaluation of Tl.nl can be found in Appendix. 



IV. RESULTS AND DISCUSSION 

In this Section, we discuss the numerical results for the leading-twist light-cone pion and 
kaon DAs with the gauged nonlocal effective chiral action. We first fix the scale parameter 
of the model, i.e. A in Eq. (|22|). In Ref. [l4T |. the A was chosen to be 1.2 GeV by using 
the normalization condition in Eq. (|19|). with which the pion and kaon decay constants are 
simultaneously reproduced. However, since the Ward-Takahashi identity for the axial- vector 
current is broken, the improper Pagels-Stokar (PS) expression for the pion decay constant 
had to be used to normalize the DAs properly. In the present gauge-invariant approach, we 
are able to determine the scale parameter A without any trouble. 

In Table |U we list the pion and kaon decay constants evaluated by using the normalization 
condition in Eq. (fTTJj) . As shown in Table HJ when we set A ~ 1.0 GeV, the values of 
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F n are in accord with the empirical one (F n = 93 MeV) for powers of the form factor 
n = 1,2,3. We find that without the nonlocal contribution from gauge invariance, the 
results are underestimated by about 10 ~ 20%. The kaon decay constant F% is found to be 
about 10 % smaller than the empirical data Fx = 113 MeV. However, it is well known that 
the meson-loop corrections play a crucial role in describing the F^ 33]. 





Ft, (Exp: 93 MeV) 


F K (Exp: 113 MeV) 


n 


Local Nonlocal Total 


Local Nonlocal Total 


1 


87.03 28.75 91.66 


92.59 31.94 97.95 


2 


84.47 30.38 89.77 


90.87 34.37 97.15 


3 


83.68 31.02 89.24 


90.31 35.31 96.97 



TABLE I: The results of the pion and kaon decay constants F n and F% [MeV]. 



In Fig. Q we depict the pion and kaon DAs in the upper and lower panels, respectively. 
The dashed, long dashed, and solid curves represent the local, nonlocal, and total contri- 
butions to the DAs, respectively. For comparison, we also draw the asymptotic DA in the 
dotted curve in each panel. As the power of the form factor n increases, the local contribu- 
tions are almost independent of n, whereas the nonlocal ones are drastically changed: Their 
humped shape is getting manifest as n increases. As a result, the total pion DA is getting 
flattened in the region of 0.25 < u < 0.75, as n increases. Although we do not present 
explicitly the pion DA with higher n, it is even more flattened, if we increase n higher than 
three. Note that the pion DA is suppressed at the end points of u. In the lower panels of 
Fig. 01 we find that while the local contribution to the kaon DA is almost symmetric, the 
nonlocal one turns out to be rather asymmetric. However, when these two contributions are 
added, the kaon DA is unexpectedly similar to the asymptotic one. However, the kaon DA 
is quite much suppressed at the end points, in particular, with higher n, compared to the 
asymptotic one. It indicates that the effects of the current quark mass becomes obvious for 
the kaon DA at the end points, as shown in Ref. The suppression of the pion and kaon 



DAs at the end points corresponds to the negative second Gegenbauer moments. 

It is also of great use to examine the Gegenbauer polynomial to analyze the DAs in 
detail 00- The Gegenbauer moments for the DAs are defined by the coefficients of the 



Gegenbauer polynomial expansion of the DAs, a m : 

(f>{u) = 6u{l - u)[l + ai Cf 2 (0 + a 2 Cl'\0 + •••]■ (25) 

Note that the odd (2m + 1) Gegenbauer moments are all zero for the pion DA because of 
isospin symmetry. The Gegenbauer moments to the sixth order (m = 6) are evaluated and 
are listed in Table |H] for the pion and in Table IIHI for the kaon. 



Our results of a\ for n = 1 are rather consistent with those of various model calculations 
listed in Table ITT1 We find that the absolute values of a£ are in general larger than those given 
in Ref. [l4[ in which the nonlocal contributions were not considered. Moreover, we verified 
that the difference between those with and without the nonlocal ones becomes evident as 
n increases. We obtain the negative values of a\ for n = 2 and 3, which are different from 
those of other models 17, 3?|. As mentioned previously, the sign for a 4 plays an important 



role in the end-point behavior of the pion and kaon DAs. The negative sign indicates that 
the DAs are suppressed at the end points, i.e. the DAs turn out to be in concave shapes at 
the end points. 
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FIG. 1: Pion (upper panels) and kaon (lower panels) distribution amplitudes. The dashed, long 
dashed, and solid curves represent the local, nonlocal, and total contributions, respectively. The 
asymptotic DA is also drawn for comparison. 



In Table ITTT1 we list the Gegenbauer moments for the kaon DA, including both odd and 
even ones. The most interesting feature shown here is that the signs of the calculated of are 
all negative. The results of the QCD sum rule calculations are not conclusive in determi ning 
the sign of af (HI 3(| 3^]. An interesting discussion on this issue can be found in Ref. |34j ^ 



In addition, the a% are also negative for n — 2, 3, being obviously different from the 
Note that the values of af are all positive for the QCDSR calculations in Refs. 3^, 36lE 



ion. 



n 


a\ 


a\ 


»6 


1 


0.11911 


0.01438 


-0.00009 


2 


0.07127 


-0.03605 


-0.02103 


3 


0.05339 


-0.06088 


-0.02596 


1 


0.56 






[17] (2.4 GeV) 




0.044 ±0.016 0.023 ±0.010 


[19] (2.4 GeV) 


0.12 ±0.03 






[34] (1.0 GeV) 


26+ - 21 

u - ZD -0.09 






[£5j (2.24 GeV) 


0.236(82) 






[36] (1.0 GeV) 


0.25 ±0.15 






[37] (1.0 GeV) 


0.44 


0.25 





TABLE II: The Gegenbauer moments for the pion DA. The numbers in the parentheses stand for 
the renormalization scale. 



S 



The analysis of the Gegenbauer moments and a\ has been performed first by Schmed- 
ding and Yakovlev [19], based on the CLEO data of the transition form factor F 7 . 7vr (Q 2 ) 
at the renormalization point /i = 2.4 GeV (SY-scale). In order to compare them to the 
present results of the Gegenbauer moments, it is necessary to evolve the scale of the present 
calculation to the SY-scale by using the one-loop QCD evolution equation. Note that the 
anomalous dimension for the Gegenbauer moments for the leading-twist DAs are obtained 
as follows 0,11 3||: 



3 + 



where 7^ = 7.11, 7^ = 11.11, 73^ = 13.95, 7 4 u; = 16.18. Thus, the Gegenbauer moments 
in two different renormalization scales, Ai = 2.4 GeV and A2 = 1.0 GeV, can be compared 
by the following relation: 



[m + l)(m + 2) 



m+l 1 

-4E 



m'=l 



(26) 



i(Ai 



t (A 2 



a(Ai 
a(A 2 



7i 0) /(2/3o) 



~ a, 



t (A 2 



ln[A 2 /A QCD ] 
lntAi/AqcD] 



(27) 



Here, we use /3 = 9 and A QCD ~ 0.2 GeV. The contours on the Gegenbauer parameter plane 
(a 2 , 04) analyzed from the CLEO experiment are depicted together with our results in Fig. 
The solid and dashed ellipses indicate the 68% and 95% confidential levels, respectively |19j . 
We also indicate the positions for the asymptotic DA (AS) and for the Chenyak-Zhitnitsky 



(CZ) DA |40j. The closed circles (•) denote the values of (a£,aj) calculated in the present 
work. On the contrary, the open circles (o) are taken from Ref. 0] i n which the gauge 
invariance of the effective chiral action was not considered. We observe that the closed 
circles for n = 1, 2, 3 are all inside the 95% confidential levels. However, when the nonlocal 
contributions are turned off (open circles), which is identical to the results of Ref. [l4|, a\ 
and a\ for n = 3 lie outside the 95% confidential regions. We also see that the values of \dl\ 
become slightly larger with the nonlocal contributions. We infer from it that the pion DA 
is getting more flattened by larger values of \a\\. 

In addition, we depict the results for A = 1.2 GeV and 1.5 GeV in the stars (*) and 
asterisks (*), respectively in Fig. |2l to test the dependence on the scale parameter of the 
model A. We see that as A increases the positions of (aj, a J) move toward the center of the 
95% confidential level. The a\ is less sensitive to the A, compared to a\. We verified that 



n 


a \ 


of 


4 


of 


n K 
a 5 


"6 


1 


-0.01879 


0.02525 


0.00467 


0.00676 


-0.00283 


0.00069 


2 


-0.01193 


-0.01935 


0.00599 


-0.02470 


-0.00371 


-0.01183 


3 


-0.01012 


-0.03694 


0.00706 


-0.03840 


-0.00566 


-0.01365 


[34J (1.0 GeV) 


07+ - 02 

u ' u ' -0.03 


27+ - 37 
u - z ' -0.12 










[36] (1.0 GeV) 




0.30 ±0.15 










[38] (1.0 GeV) 


-0.18 ±0.09 0.16 ± 1.10 










[22] (1.0 GeV) 


0.10 ±0.12 













TABLE III: The Gegenbauer moments for the kaon DA. The numbers in the parentheses stand for 
the renormalization scale. 
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Gcgcnbaucr parameter space (fx = 2.4 GeV) 




FIG. 2: Gegenbauer parameter space (02 — 04) for the pion DA with the CLEO experimental data 
analysis taken from Ref. [IsL in the renormalization scale /x = 2.4GeV. The solid and dashed 
ellipses indicate the 68% (la) and 95% (2a) confidential levels of the analysis. The open circles 
(o) denote the results of the Gegenbauer moments without the nonlocal contributions (A = 1.2 
GeV) 14]. The close circles (•), stars (★), and asterisks (*) stand for the cases with the nonlocal 
contributions for A = 1.0 GeV, 1.2 GeV and 1.5 GeV, respectively. The results with the same 
power n and different A are connected by the solid line. 



the calculated pion decay constants F n for A = 1.2 GeV and 1.5 GeV are only smaller by 
about 10% than the empirical one. 

To show the dependence of the pion DAs on the scale parameter A, we draw in Fig. El the 
pion DA with different A. We observe that the pion DAs are more suppressed at the end 
points and more humped as A increases. 



A = 1.0 GeV to jt = 2.4 GeV A = 1.2 GeV to fi = 2.4 GeV A - 1.5 GeV to fi — 2.4 GeV 




FIG. 3: Pion distribution amplitudes at fj, = 2.4 GeV with A = 1.0 GeV, 1.2 GeV and 1.5 GeV 
from the left end panel. 
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V. SUMMARY AND CONCLUSION 



We investigated the pion and kaon distribution amplitudes within the framework of the 
nonlocal chiral quark model from the instanton vacuum, making the gauge invariance of 
the low-energy partition function preserved in the presence of the external fields. The 
scale parameter of the present work was set to be A ~ 1.0 GeV which is not far from 
that determined by using the saddle-point equation. Using A ~ 1.0 GeV, we were able to 
reproduce the empirical values of the pion and kaon decay constants qualitatively well by 
using the normalization condition of the distribution amplitudes. The effects of flavor SU(3) 
symmetry breaking effects were explicitly taken into account consistently. 

We found that the nonlocal contribution arising from the gauge invariance of the effec- 
tive chiral action makes the pion distribution amplitude flattened more in the region of 
0.25 < u < 0.75, compared to the former investigation without consideration of the gauge 
invariance |l4j . The results indicate that the absolute values of Gegenbauer moments a.2 and 
04 become larger. As for the kaon, the nonlocal contribution is less important. We noticed 
that af is still ne gati ve even with the nonlocal contributions, which is similar to those in 
the previous work ^ . 



By virtue of the renormalization group equation, we evolved the Gegenbauer moments 
ci2 and CI4 so that we may compare them with the Schmedding-Yakovlev analysis on the 
CLEO experimental data at fi = 2.4 GeV. We found that the present results with the gauge 
invariance turn out to satisfy the Schmedding-Yakovlev analysis of the Gegenbauer moments 
at the 95% confidential level. We found that the values of (af , af ) were very close to those 
of the asymptotic distribution amplitude. 

In conclusion, the present gauge-invariant approach has several advantages, compared 
to the former works. Firstly, it provides a more reasonable scale to the pion and kaon 
DAs, since we can use the normalization condition without any gauge-invariant problem. 
Secondly, the nonlocal corrections arising from the gauge invariance of the effective chiral 
action play a significant role in describing the pion distribution amplitude. Compared to 
the Schmedding-Yakovlev analysis, the present results are phenomenologically better than 
those of the former work. In order to describe the pion and kaon distribution amplitudes 
more consistently, we want to introduce the meson- loop corrections {1/N C corrections) to 
them. The corresponding work is under progress. 
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Appendix 

The first term in the square parenthesis in Eq. (|24j) reads 



tr„ 



M f (k) ^M f (k-P) 

fij5 t^TZ 7^75 



D(k) 



D(k-P) 



-AM Jf(m f )f(m f ) 



nA 2 



nA 2 — k 2 



nA 2 



nA 2 - (k — P) 2 
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x <k ■ n 



m g + M f(m g 



nA 2 



nA 2 — (k — P)' 



2n 



(k-P)-n 



m f + M f(m f ) 



nA 2 


2n 




nA 2 — k 2 




! 



X 



k 2 — m 2 — 2M mff(mf) 



nA 2 


2n 


-M 2 f 2 


nA 2 


4n^ 


nA 2 — k 2 




nA 2 — k 2 





x \(k-P) 2 -m 2 -2M m g f(m g ) 



nA 2 



2 ii 



- M 2 f 2 



nA 2 



nA 2 - (k - Pf 



x 



nA 2 - (k-P) 2 _ 
4V^(o*_ - 7 /) B (/%- - 7,)" 

a(ak_ - lf ) 2n [m g ((3k_ - lg ) 2n + Vg ] - (3(f3k_ - lg ) 2n [m f (ak_ - lf ) 2n + Vf ] 
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where 



i.e. 



V f = (ak- - D f )(ak- - 7/) 4n - 2r) f m f (ak- - 7/) 2n -r]j + 
V g = (Bk- - D g )(Bk- - lg f n - 2 Vg m g (Bk- - lg ) 2n - r, 2 + ie. 

We use the following parameterizations for convenience: 

(k - P) 2 = (u- l)P+k- -(u- l)ml -k 2 T = -{u- \)m\ - k 



Ti 



nA 2 - k 2 = -uP + k- + [k^, + nA 2 ] = -ak_ + 7/, 

nA 2 - {k — P) 2 = -{u - l)P+k- + [{u - \)m\ + k 2 T + nA 2 } = -/3k- + 7 9 , 
k 2 — m 2 r = ak- — [k^ + m 2 A = ak- — 5 



/• 



(k - P) 2 -m 2 g = (3k_ - [(u - l)mj + k 2 T + m 2 ] = /3Jfe_ - 5. 
M f(m f )(nA 2 ) 2n = Vf , M f(m g )(nA 2 ) 2n 
M f (k) - Vf 

1 



'91 



(ak- — 7/) 2ti 

(ak- - 7/) 4 " 

D(k) V f 

1 (Ph- - l 9 ) 



M f (k - P) 

+ 771/ + 



Vy 



((3k_- lg f^ 
Vf 



-In 



D(k-P) 



V n 



(ak-.-jf) 



2n 



Vf 



((3k- - lg ) 



2n 



An 



The second term can be evaluated using the dynamical quark mass of Eq. 
2ntr^ 



n 



kJMf(k)JM f (k-P) 



D(k)(nA 2 - k 2 ) 
x \m f + M f(m f ) 



nk ■ hM J f(m f )f(m g 



(nA 2 - k 2 ) 

k 2 -m 2 f - 2M f(m f ) 
8na y /rjjrj g ~ [mf(ak- — 7/) 2n + Ej\ (ak- — 7/)™" 1 



nA 2 



nA 2 - k 2 



as follows: 

nA 2 



nA 2 - (k - Pf 



nA 2 


2n^ 


nA 2 — k 2 





nA 2 


2n 


- M 2 f 2 


nA 2 




nA 2 — k 2 




nA 2 — k 2 





-1 
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Similarly, the third term can be written as follows: 

^iPy/Wb [rrigiPk- ~ lg? n + %\ (£k- - i g ) 
" " {uk-- lf ) n V g 

For convenience, we denote these three terms evaluated above by and T^ L , which is the 
sum of the second and third terms. 



2ntr^ 



(k-P)-fiy/M f (k)y/M f (k-P) 

D(k-P)[nA 2 - (k-P) 2 ] 
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